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Temperature dependence of the thermoelectric power (TEP) of metallic systems with cerium and 
ytterbium ions exhibits some characteristic features, which can be used to classify these systems into 
distinct categories. We explain the observed properties in terms of the Kondo effect modified by 
the crystal field (CF) splitting and relate various shapes to different energy scales that characterize 
scattering of conduction electrons on Ce and Yb ions at different temperatures. The low- and 
high-temperature behaviors are calculated for different fixed point models and the overall shape of 
the TEP is obtained by interpolation. At high temperatures we use the Coqblin-Schrieffer model 
(CSM) and calculate the thermoelectric power by the perturbation expansion with renormalized 
coupling constants. Performing the renormalization by the 'poor man's scaling' (PMS), we find for 
large CF splitting two Kondo scales, Tk and 3> Tk- The thermopower obtained in such a way 
exhibits a large maximum at about , and a sign change at about Tk- The PMS also shows that 
for T <,T§ the / level behaves as an effective multiplet with the degeneracy of the CF ground 
state and Kondo scale Tk- Thus, we assume that the low-temperature fixed point models do not 
require the CF splitting. The properties of dilute Ce and Yb alloys are obtained for T < Tk from 
an effective spin-degenerate single-impurity Anderson model (SIAM). The parameters of SIAM are 
such that its effective Kondo scale coincides with the Kondo scale Tk of the CSM. The stoichiometric 
compounds are described by an effective spin-degenerate periodic Anderson model (PAM), such that 
its characteristic energy scale To is the same as Tk- The transport coefficients of PAM and SIAM 
are calculated by perturbation theory. The interpolation between the Anderson model results, valid 
for T < Tk, and the CSM, results valid for T > Tk, provides a qualitative explanation of the TEP 
of most Ce and Yb intermetallics. 
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Introduction 

The thermo-electric properties of intermetallic com- 
pounds with Ce and Yb ions with one / electron or / 
hole show many interesting features and their thermo- 
electric power (TEP) considered as a functios, of temper- 
ature, S{T), assumes a variety of shapesB~E3 The func- 
tion S{T) is often non-monotonous and in some systems 
it exhibits a sign change at lowest temperatures. Above 
100 K, the TEP can assume giant values, and much of 
the recent interest in heavy fermion thermo-electricity 
is due to the belief that some of the new systems, with 
the thermopower. larger than 150 fiV/K, might be useful 
for application.tj The functional form of S{T) appears 
to be quite complicated, but systems with similar ther- 
mopower exhibit similarities in other thermodynamic and 
transport properties as well. The shape of S{T) can be 
used to classify jCe and Yb intermetallics and alloys into 
distinct groups .l3 

The thermo-electric anomalies due to the care-jsarth 
ions have stimulated a lot of theoretical workEJcj It is 
now clear that most of the experimental results can be 
understood in terms of the exchange scattering of conduc- 
tion electrons on the 4f^ state of Ce or the 4f^'^ state of 
Yb provided one takes into account the splitting of the / 
states due to the crystalline electric field (CF) . The com- 



plete description would have to consider the many-body 
effects due to the local correlation and the hybridization 
of / electrons with the conduction states, together with 
the point group symmetry of the crystal and the spin- 
orbit and the CF splittings. Although the solution of 
such a general model is not available, one can assume, us- 
ing the usual renormalization group arguments, that the 
low- and high-temperature behaviors belong to different 
fixed points. Provided we know the solution of various 
models pertaining to different fixed points we can then 
construct the overall temperature dependence of S{T) by 
interpolation. Here we show that the typical shapes of 
the thermopower found in metallic systems with Ce and 
Yb ions can be explained in such a simple way. The 
magnetic response of these systems has been explained 
recently using the same approach. 

At high temperatures we use the Coqblin-Schrieffer 
model (CSM), which takes into account the degener- 
acy and the CF splitting of Ce 4f^ electron and Yb 
4f^'^ hole. Thus, we describe the / ions as local mag- 
netic moments and assume that the scattering of con- 
duction electrons on such moments is incoherent regard- 
less of their concentration. This is consistent with the 
observation that the high-tcmpcrature thermopowcrs of 
dilute alloys and atpichiometric compounds are often sur- 
prisingly similar .113 The giant values of S{T) around the 



1 



high-temperature maximum are explained very well by 
the 3rd-order rpeiiuijbation theory for the transport re- 
laxation time.E3ci'E§E3 The perturbation theory for the 
CSM can be further improved by renormgdiziag the cou- 
pling constants-by the poor man's scalingElO The scal- 
ing solutionE3~a shows that a large CF splitting A leads 
to two Kondo temperatures, Tk and ^ Tk- The 
high-temperature maximum of S{T) relates to Tf and 
the reduction of temperature can lead to a sign change 
at Tk- Typical values of Tk,T^, and A in Ce and Yb 
system are 10 K, 100 K, and 100 - 1000 K, respectively. 
However, the renormalized perturbation expansion be- 
comes unreliable around Tk and cannot be used to de- 
scribe the non-monotonic thermopower found in many 
heavy fermions below Tk- 

To obtain the fixed point models appropriate for the 
description of the system at low-temperatures, we use the 
scaling result that the CF split / state of the CSM can be 
approximated for T <C T^ by an effective multiplet with 
the degeneracy of the CF ground state. Thus, as regards 
the low-energy dynamics, we assume that the sole effect 
of the excited CF states is to provide the right Kondo 
scale for an effective model. Once the lowest energy scale 
is set to Tk, the excited / states can be removed from the 
low-temperature problem. In addition, we assume that 
for T ~ Tk the scattering of conduction electrons is inco- 
herent, regardless of the concentration of the rare earth 
ions. The properties of the system for T < Tk arc then 
found by considering the 4f ions as an assembly of /levels 
with the degeneracy of the CF ground state and no CF 
splitting. The transport and the thermodynamic anoma- 
lies are related to the scattering of conduction electrons 
on such effective / states. However, to allow for the possi- 
bility that for large concentration of Ce and Yb ions the 
ground state of the system could become coherent, we 
treat the dilute alloys and the stoichiometric compounds 
in a different way. 

In dilute alloys, where the scattering of conduction 
electrons on rare earth ions remains incoherent down to 
lowest temperatures, the effective low-temperature prob- 
lem for the Ce and Yb impurities is described by the 
single-impurity Anderson model (SIAM). The effective 
flevel of the SIAM is assumed to have the same degen- 
eracy as the CF ground state of the 4/^ level and the 
effective parameters are adjusted in such a waff that the 
width of the it f-electron Kondo resonanceEil coincides 
with the Kondo temperature Tk of the CSM. Thus, we 
assume that for T < Tk the CSM can be represented by 
the SIAM. The TKP-pf the SIAM is then found by some 
accurate methodoa 

In stoichiometric compounds, the problem is somewhat 
more difficult, because the reduction of temperature leads 
to the crossover to the coherent regime. We assume that 
for T < Tk the properties of the lattice of Ce or Yb ions 
are described by a periodic Anderson model (PAM) , such 
that the effective / levels have the same degeneracy as the 
lowest CF state of the 4f^ or 4f^'^ electrons. The absence 
of the excited CF states simplifies the calculations and 



the low-temperature transport coefficients of the PAM 
can be obtained by the self-consistent perturbation the- 
ory with respect to the correlation UxB This leads to 
a characteristic energy scale Tq, such that for T < Tq 
the model describes a coherent Fermi-liquid state, and 
for r > To an incoherent state. That is, above Tq, we 
consider the rare earth ions as independent scattering 
centers for conduction electrons. The effective parame- 
ters are now adjusted to give Tq ~ Tk, and we assume 
that for T <Tk the PAM and the CSM are dynamically 
equivalent. Above Tk, the physical state begins to de- 
viate from the incoherent high-temperature state of the 
PAM which neglects the excited CF states. For T > Tr: 
the low-energy properties of stoichiometric compounds 
are described by the CSM with the CF splitting. Once 
we obtain the solution of various fixed point model, the 
overal temperature dependence is found by interpolation. 

The paper is organized as follows. First we describe 
the experimental results. Then we present the scaling 
solution of the CSM and discuss the TEP obtained by 
the renormalized perturbation theory. Next, we summa- 
rize the low-temperature results obtained for the SIAM 
and the PAM. Finally, we obtain the full TEP by inter- 
polating between various fixed point solutions and use 
these results to discuss the data. 



Classification of the experimental data 

The experimental results for TEP of the Ce- and Yb- 
based intermetallic compounds exhibits some character- 
istic features, which can be |Used to divide these systems 
into several distinct groups.t3 

The type (a) thermopower is characterized by a deep 
negative minimum at low temperatures and a broad pos- 
itive high-temperature peak between 100 K and 200 K. 
The type (a) behavior is found in a large majority 
Ce-ba^e^ heavy fermi|Ojj-C|pajipounds like.jCjeCu2Si2,l 

CePd2Si2,ElE3 CeP|^n,[ 



CePb3llII^CeCu2 
CePdGe,EiCeAl 



17,25U28 



CeRh2-zNi2Si2 for small z and 



CePtGe2.i--J It is also found in Yb-based systems like 
YbAgCu4 or YbPd2Si2El 

The type (b) thermopower is similar to the type 
(a), except that at very low temperatures the ther- 
mopower changes sign again and exhibits an addi- 
tional positive peak. These features, are often fompd, 
in dilute Ce-allovs, like Cea;Lai_a;Ni|13 Cea;Lai_a;Al3,cll 
Cea;Lai_a:Pd2Sj2p3 and in some concentrated systems 
like Cea:Lai_a;a and Ce_Lai_2;Cu2Sz2j3 or jStoichiomet- 
ric compounds CeAlsJsllJ or YbAgCu4.Ell The ad- 
ditional low- temperature peak emerges in CePd2Si2 
with pressjure,Ea and in Cea;Lai_2;Pd2Si2 with chemical 
pressure.ES In some cases the positive low-temperature 
peak is concealed by the supercondueting or magnetic 
transitions. For example, in CeCu2Si2El the positive low- 
temperature upturn does not appear in the zero-field 
data and it is only seen in an external magnetic field 
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which suppresses the superconducting transition. 

In type (c) systems the low-temperature maxi- 
mum is more pronounced than in type (b) but 
the separation to the high-temperature maximum is 
smah and such that the minimum in-between does 
not extend down to negative values. The sign of 
S{T) remains in type (c) systems the same at all 
temperatures. These features ara found in dilute 
systems like Cea;(LawY2)i_2;Al2,y Ce^Lai_a;Cu6,t3 

and Ccj^Lai- ^U2Si2,llil ilH 

Ce-rich compounds like 

for smaJl 



^Pbi_^Sn^),El CcRu2Si2,EI Ce(Cu^Aui_ 



zfH and 

Ce2,Lai_a;Cu2Si2 



-a:Nio.8Pto.2,E3 



j in Ce 
I Ce2:Yi_2:Cu2Si2 
The type (d) thermopower exhibits a large high- 
temperature peak, and perhaps a shoulder on the low- 
temperature side of that peak. This behavior is found 
at ambirait pressure in CeCug and GexLai-xCug, — 
CeInCu2,El CeCuaGaa, and CeCugAlaM CePNiGcE 
and under Jiigh pressure in CeCu2Si2,El CeCu2Ge2,c3 and 
CePd2Si2.El 

Finally, the TEP of shape (e) is found in valence fluc- 
tuators, where one observes a monotonic increase all t 
way up|-tp room teroaeratures. This is seen iruCeBeia 
CePdaP CeNi25'i2,HB or in CeNi and CeNi2.Ll 

The border between various categories is not very 
sharp and the TEP can be transformed from one type 
to another by changing the environment of 4f ions. 
In addition to thet-already listed examples, we men- 
tion Ce2;Lai_2;Cu6,c3 where for x < 0.5 one finds two 
peaks separated by a well resolved minimum, while for 
X > 0.9 only a sinsle hump with a shoulder on the 
low-temperature sidd£j remains. Similarly, by chang- 
ing the concentratkm of Ce ions, one transforms S{T) 
of Ce2:Lai_a;Pd2S|iaE3-i'mm type (a) to type (b), that of 
Ce2,Lai_:rCu2Si2anJ't3'K3,from type (a) to type (c), and 
that of Cea;Yi_:rCu2Si2E3 from type (a) to type (d). Di- 
rect application of pressure (rather thap-chemical pses- 
sure) ti^sforms the TEP of CePd2Si2,Ea CeCu2 812^1 or 
Ce AI3JI3 from type (a) or (b) to type (c) or (d). 

Typical behaviors of the TEP are found in the 
Cei_a;Laa;Cu2Si2 family of intepnetallic compounds for 
0.01 <C a; < 1.0 at % of Ce^ where the TEP data 
exhibit nearly all the shapes discussed above. (The pos- 
itive low-temperature peak which is seen in the dilute 
limit is also found for large concentration of Ce ions, 
provided one suppresses the superconducting transition 
by the magnetic field.) The peaks in the thermoelec- 
tric power correlate with the logarithmic behavior in tha 
electrical resistance; the magnetic susceptibility dataO 
indicate that above 100 K the local moment is 6-fold de- 
generate, while below 20 K it appears to be only 2-fold 
degenerate. 

In summary, the experimental data show that the char- 
acteristic shapes (a) - (e) appear regardless of the con- 
centration of magnetic ions. The high-temperature ther- 
mopower is insensitive to small changes in the concentra- 
tion of rare earth ions and the thermoelectric anomalies 
are accompanied by the anomalies in other transport and 



thermodynamic properties. We take this behavior as an 
indication that the TEP anomalies are due to the scat- 
tering of conduction electrons from the 4-f state of Ce 
or Yb ions, and that the high-temperature scattering is 
incoherent regardless of the concentration of magnetic 



High-temperature approximations 

The high-temperature properties of the rare earth ions 
in metallic environment are modeled Coqblin-Schrieffer 
HamiltonianeS with the CF splitting. 
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where all the symbols have their usual meaning. The first 
term describes the CF-split 4f^ state of Ce ions or 4f^^ 
state of Yb ions, the second term describes the conduc- 
tion band of width 2Z?o and a constant density of states 
po, the third term defines the pure (non-diagonal) ex- 
change scattering between 4f^ states and band electrons, 
and the last term is the total (diagonal) potential scat- 
tering. We consider here the antiferromagnetic coupling 
only. The summation over v is over all the CF states, 
and Jo and I^) are the initial coupling constants, which 
are assumed to be i^- independent. For simplicity, we rep- 
resent the 4f^ electrons of Ce (the 4f^'^ hole of Yb) by 
their lowest J = 5/2 (J = 7/2) spin-orbit state, and con- 
sider the CF scheme in which the ground state level at 
Em and the excited state at Em are a,„-fold and 
fold degenerate, respectively. The energy separation is 
Em — Em = A. The model neglects the coherent scatter- 
ing on Ce ions and is most appropriate for dilute alloys. 
However, it also applies to concentrated systems at tem- 
peratures such that the scattering of conduction electrons 
on the rare earth ions is incoherent. 



Perturbation theory for the thermoelectric power 

We describe the high-temperature heat and charge 
transport by the Boltzmann equation and evaluate the 
scattering rate for |tbe-.tcansport relaxation time up to 
the 3rd order in Jo.CJ'l3o All the computational details, 
can be found in the papers bypCornut and CoqblinE^I 
and Bhattacharjee and Coqblinj23 which we refer to as 
CC and BC, respectively. Here we just quote the ther- 
mopower result obtained by BC for the doublet-quartet 
CF scheme {am = 2, aM = 4), 
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where Sa and i?A are given by 
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The derivative of the psi-(digamma-) functionEl is de- 
noted by ip'{x), the occupancy of the CF states is 
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and the effective potential scattering is 

Vm = Vm- Jm {rim) , Vm = Vm ~ Jm ("a/) 



these features follow from the perturbation expansion, 
provided we renormalize the coupling constants by the 
poor man's scaling. 

The scaling approach to the single-impjyititji, models is 
j^g-j explained in great detail in many paper^To and here 
we just use these results to calculate the thermopower. 
The scaling equations are generated by reducing the con- 
duction electron cutoff from Dq to D and simultaneously 
rescaling the coupling constant, J{D), so as to keep tke. 
low-energy excitations of the total system unchanged^ 



(4) 



(5) 
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Note, poSa/Ra is a dimensionless quantity, po is the den- 
sity of conduction states arround the chemical potential 
/i, and the prefactor in Eq. (||) is ks/e « 86ijlV/K). 

The qualitative features of the thermopower described 
by Eq. (|) with constant and isotropic coupling constants 

(Jm = Jm = JraM = Jq, and Vm = Vm = V) follow 
straightforwardly from the asymptotic expansion of the 
elementary functions in Eq. (Q) and can be summarized 
as follows (detailed analysis can be found in Ref. BC). 
At temperatures such that T > A the thermopower is 
small and behaves as l/T. The sign of S and its slope 
depend on the relative size of Jq and V. As tempera- 
ture decreases, the thermopower increases up to a large 
(positive or negative) value at about Tmax — A/3. At 
temperatures such that T <C A we find Gi oc T/A and 
S{T) cx T. 

For parameters in the physical range one easily finds 
the peak value of thermopower above 50 /iV/ K and 
Eq. (0) captures the essential high-temperature features 
of Ce and Yb intermetallics.Ej At low temperatures, how- 
ever, perturbative approach leading to Eq. (|^) fails and 
the observed, non-monotonic low-temperature behavior 
is not reproduced. 



Poor man's s 



The change of sign of the TEP below Tmax, which 
is often seen in the experimental data, cannot be ob- 
tained from Eq. (H) for any value of Jq and Vq. However, 



The solution for the CSM, va 
renormalizcd couplings, readsl 



up to the 2nd order in 
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where Tx is the Kondo temperature, defined by Jq and 
Dq. For a given set of initial parameters the renormalized 
coupling constant J{D) is completely characterized by 
Tk- The result obtained for m = 2, M = 4, and A = 330 
K is shown in Fig. |^ as a function of D/T^, where T^ is 
defined below (see Eq. |l^) . For a given value of Dq and 
A, the Kondo temperature grows exponentially with Jq; 
for a given Dq and Jo the Kondo temperature grows with 
the CF splitting as (A/Dq)^'^. The scaling trajectory 
defined by Eq.(||) diverges at D = Tk, such that J{T) 
can be used to renormalize the weak-coupling theories 
for D ^ Tk only. Note that the CF splitting lowers the 
degeneracy of the / level and reduces the Kondo scale 
from T^ to Tk, which extends the validity of the weak- 
coupling theory to much lower temperatures. 
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FIG. 1. D dependence of the scaled exchange coupling 
J{D) of CS Hamiltonian in the case of a ground state dou- 
blet and an excited quartet, separated by the CF splitting A. 
Solid curve is the 2nd-order Hanzawa scaled result. Dashed 
and dot-dashed curves are the high-energy and low-energy 
asymptotic scaling trajectories appropriate to a sixfold de- 
generate and a doubly degenerate level, respectively. 

The scaling trajectory J{T), defined by^^q. ^) and 
shown in Fig. |l|, has two asymptotic regimes.EJ The higli- 
temperature asymptote, Jh{T), is obtained from Eq. (|^) 
by neglecting A with respect to D, which leads to 
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where is the high-temperature Kondo scale given by 



(Tf)^ = (TKr(rK + A: 



A// 



(10) 



and N = m + M. The same trajectory would be obtained 
for an effective sextet without the CF splitting, assuming 
the initial condition Jh{Dq) — Jq. 

The low-temperature asymptote Jl(D) is obtained by 
neglecting in Eq. Tk and D with respect to A, i.e., 
by neglecting the effects of the excited CF levels. This 
identifies Jl{D) as the scaling trajectory of an effective 
doublet with Kondo scale Tk and the effective coupling 
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Note that the value of Tk in the above expression is 
defined by Eq. (^) with the initial condition J (Dp) — Jq. 
The initial condition Jl{Do) ~ Jo smd Eq. ( |l l| ) would 
lead to a Kondo temperature T^ that would be much 
lower than Tk, i.e., the Kondo scale that enters Eq.([ll|) 
is exponentially enhanced with respect to Tj^ because of 
CF splitting. 

Thus, the scaling analysis shows that the / level be- 
haves at high temperatures as an effective sextet with 
the scaling trajectory Jh{D) and Kondo scale T^ , while 
at low temperatures it behaves as an effective doublet 
with the scaling trajectory Jl{D) and Kondo scale = 
Tk- The scaling trajectory J{D) of the CSM with the 
CF splitting interpolates between these two asymptotic 
regimes. Here, we discussed the scaling solution obtained 
by the 2nd-order perturbation expansion butj-the same 
behavior is obtained in the next leading order.H 



Renormalized perturbation expansion 

The scaling theory can be used to obtain a simple de- 
scription of the overall magnetic and transport properties 
of the CSM. It provides the thermopower in the following 
way. We reduce the conduction-electron half-bandwidth 
from Dq to D — AksT, where A is a numerical con- 
stant of the order of unity, and find the effective couplings 
J(T), Vm{T), and Vm{T). Next, we notice that the re- 
duction of the bandwidth does not change the form of the 
Hamiltonian, and that the form of the response functions 
obtained by perturbation expansion in terms of the effec- 
tive coupling constants remain invariant with respect to 
the scaling. Thus, to obtain the poor man's scaling re- 
sult for the thermopower we introduce the temperature- 
dependent cutoff A ksT, find J(T) from Eq.®, and 
substitute this renormalized expansion parameter in the 
expressions (§) - (0). 

To apply the scaling result to a specific Ce or Yb sys- 
tem we have to determine the degeneracy and the split- 
ting of the CF levels, estimate the initial coupling con- 
stant (or Tk), and choose the cutoff constant A. The 



thermopower obtained by the renormalized perturbation 
theory depends on the magnitudes of the Kondo tem- 
perature, the potential scattering, the CF splitting, and 
the density of states of the conduction electrons at the 
Fermi level. These parameters are restricted to a rather 
narrow range by the requirements that the model explain 
not just the thermopower data but that it also leads to a 
consistent description of other thermodynamic and trans- 
port data of a particular system. We estimate the CF 
splitting and the degeneracies of the CF states by ana- 
lyzing the neutron scattering and the high-temperature 
magnetic anisotropy data. The Kondo temperature Tk, 
i.e., the initial coupling Jq, can be obtained from the 
magnetic susceptibility data, provided we know how to 
extract^he single-ion contribution from the experimental 
resultal3. The high-temperature Kondo scale T^ is then 
found from Eq. (|lO|). The relationship between Tk and 
J(T) depends to some extent on the scaling procedure 
and the value of the cutoff constant A. Our analysis of 
the thermopoweiLpf the CSM is based on the 2nd-order 
scaling equationEHl and we j.iae, for simplicity, A = 1. The 
3rd-order scaling equations and/or some larger values of 
A do not change the TEP in any qualitative way. The 
potential scattering shows up in the electrical resistance 
and typical results for the high-temperature phaseEJ in- 
dicate rather large values for V. 
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FIG. 2. The thermopower is plotted as a function of re- 
duced temperature T /Tk and for various values of the po- 
tential scattering. Dot-dashed curve poV — —0.15, dotted 
curve poV = —0.20, dashed curve poV = —0.35, solid curve 
poV = -0.5. Here we used Tk = 8 K and A = 350 K, such 
that r# = 99 K. 

The thermoelectric power obtained for the doublet- 
quartet CF scheme, with A = 350 K and for several 
values of Tk and V, is shown in Figs. |2| - ^ The density 
of conduction states, which enters the thermopower cal- 
culations, is set to po = 2.2 (cV)~^. Fig. |^ shows the 
effect of V on S{T) for = 8 K, = 99 K, and 
A = 350 K. The details of the shape of S{T) depend in a 
rather complicated way on the relative magnitude of Tk, 
V, and A. For large V the value of T^ax is mainly deter- 
mined by T^ , which is related to Tk and A by Eq. ([l^). 
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For small values of V, the TEP changes sign before the 
full maximum can develop. The crossing temperature T^, 
where S{T) changes sign, i.e. S{Tx) — 0, increases as \V\ 
decreases. 




200 
T[K] 

FIG. 3. The thermopower is plotted as a function of tem- 
perature for poV — —0.35 and A = 350 K, and for three val- 
ues of Tk- Dotted curve = 2 K, Tf = 63 K, dashed curve 



Tk 



K, = 99 K, solid curve Tk = 32 K, 



158 K. 



Fig. I shows S{T) for A = 350 K, f = 0.35, and for 3 
values of Tk- For such a large value of po^, we find that 
the maximum of S{T) is around and that Tx scales 
approximately with Tk- Fig. ^ shows the effect of A on 
S{T) for T/f = 8 K and y = -0.35. The position of T^ax 
is rather well described by T^ , while the value of T^ is 
not much changed by A, as it mainly depends on Tk- 
The overall features of S(T) shown in Figs.^ - HfoUow 
straightforwardly from expressions (||) - (Q) . Eq. con- 
tains the products of logarithmic functions and at high 
temperatures, where J(r)pGhanges very slowly, S{T) re- 
duces to the BC expressiorO. Eq. (H) shows that the sign 
of the thermopower is determined by the temperature de- 
pendent factor J(T) — V, and that S{T) approaches zero 
for J{T) ~ V- 

However, to make a quantitative analysis of S{T) close 
to Tx, we should not neglect the higher-order terms in the 
expansion of transport integrals. The 3rd-order pertur- 
bation theory provides an indication that S(T) changes 
sign below T^ax and exhibits non-monotonic features. 
Since S(T) must vanish at T = 0, the finite value of T^ 
must lead to a low-temperature minimum, which is here 
obtained from the exchange scattering, without invok- 
ing the interaetion effects or long-range antifcrromagnetic 
fluctuations. t3 




FIG. 4. S{T) is plotted as a function of reduced temper- 
ature T /Tk for p^V = —0.35 and for various values of the 
CF splitting. The curves are obtained for Tk = 8 K, and 
the dot-dashed, solid, dashed, and dotted curve corrspond 
to A = 50K (Tf = 27 K), A = 200 K (Tf = 68 K), 
A = 350 K (rf = 99 K), and A = 500 K (Tf = 126 K), 
respectively. 

The properties of the model for T < Tk can not be ob- 
tained from scaling. On the one hand, Eq. is derived 
by assuming that thermal fluctuations do not excite too 
many conduction electrons or holes up to states near the 
(effective) band edges; this assumption does not hold for 
T < Tk, i.e., the scaling breaks down when the effec- 
tive (renormalized) band edge is of the order Tk- On 
the other hand, the perturbation theory for the response 
functions breaks down close to Tk, because the effective 
coupling constants become too large for the lowest-order 
perturbation expansion to be valid. However, Fig. |l| in- 
dicates that the breakdown of scaling occurs at temper- 
atures such that J(T) ~ Jl{T), and that the low-energy 
dynamics of the CSM is rather well described by an effec- 
tive doublet or quartet. Thus, we neglect the CF split- 
ting at low temperatures and calculate the thermopower 
of Ce and Yb intermetallics by using the methods that 
are more accurate than scaling. 



Low-temperature approximations 

To discuss the low-temperature properties we start 
from the scaling result that the 4/^ or 4/^'^ states of 
the rare earth ions behave for T ~ <C as effec- 
tive levels with the degeneracy of the CF ground state. 
We assume that the coupling of such effective states to 
the conduction electrons characterizes the Kondo tem- 
perature which is the same as Tk of the full CS model, 
and describe the low-temperature data by models with- 
out the CF splitting. Since the reduction of temperature 
could lead to coherent scattering of conduction electrons 
on rare earth ions, stoichiometric compounds are treated 
differently from dilute alloys. 
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Periodic Anderson model 

In ordered compounds, we assume that the coherence 
temperature of the lattice of rare earth ions, Tq, is not 
much different from Tr-, the lowest Kondo temperature 
of the CS model. Since the crossover from the full CF 
multiplet to an effective low-temperature CF state takes 
place much above Tr-, we neglect the CF splitting in the 
coherent regime. Thus, for T < Tq we describe the sto- 
ichiometric compounds by the periodic Anderson model 
(PAM) which has the degeneracy of the lowest CF state 
and adjust the effective parameters in such a way that 
the low-energy scale of the PAM coincides with Tk of the 
CSM. Here we use the twofold degenerate PAM defined 
by the Hamiltonian 
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(12) 



A; <T 



Cka^i<7 + h.h.) + Uy^^a\^ai^^a\^ai^, 



k ia 



where a labels two spin states of the localized level, Ej 
is the unrenormalized position of the effective CF ground 
state, V is the hybridization, and U is the Coulomb re- 
pulsion. We consider the case where the / level is below 
the Fermi level, and the width of the / band W is smaller 
than the on-site correlation. 

A simple but non-trivial approximation to study the 
PAM is provided by the 2nd-order perturbation the- 
ory (SOPT) with respect to the correlation U . The 
SOPT, which may be regarded as the simplest extension 
of Hartree-Fock theory, properly reproduces Fermi liq- 
uid properties and mass enhancement, and it reproduces 
a characteristic low-temperature scale Tq. A disadvan- 
tage of the SOPT is the fact that this Tq has not the 
proper (non-analytic) dependence on the model param- 
eters (hybridization V and correlation U, in particular), 
as it should according to the Schrieffer- Wolff transforma- 
tion. But on the other hand, arbitrarily large mass en- 
hancements and any value of Tq can be achieved in SOPT 
by choosing appropriate values for the parameters V and 
U . Therefore one can choose the parameters of the effec- 
tive PAM in such a way that To corresponds to the true 
Tk of the CS model. The SOPT is moat^easy to apply 



in the limit of high dimensions d ^ cxd,E_I because then 
the self-energy is k-independent (site-diagonal), which is 
already a good additional approximation for the realistic 
dimension d = 3. Consistent with the site diagonality 
of the self-energy is the vanishing of vertex corrections. 
Then the transport coefficients can be calculated from 



the Kubo formu 
tivity is given byl: 



The static (zero-frequency) conduc- 



2TTn 



(13) 



where a is the lattice constant, d the dimension, t the 
hopping matrix element of the band electrons, f{E) the 
Fermi function, and the function L{E) is given by 



L{E)^-Y.^lmGUE + ^0)f 



The function 



V2 



-Ef-T,(z) 



£k 



(14) 



(15) 



is the band electron Green's function of the PAM, and 
S(z) the it f-electron self-energy, which we determine ap- 
proximately within SOPT. The TEP is also determined 
by the function L{E) according to 



S = 



df_ 

dE 



{E~p.)L{E) 



eTjdE[-^)L{E) 



(16) 



As mentioned above, the SOPT yields a characteristic 
low-temperature scale Tq which can be determined or de- 
fined as follows. The it f-electron density of states (DOS) 
at the chemical potential /i is strongly temperature- 
dependent and decreases as temperature increases on 
the scale Tq towards an asymptotic T-independent value. 
Therefore, Tq can be defined as the half-width of the T- 
dependent part of the it f-DOS at fi. When calculating 
the T-dependence of the resistivity R{T) = 1/ axx from 
(p^), one obtains for most choices of the parameters the 
following characteristic behavior: A residual resistivity 
approaching zero for T — > 0, a dependence for very 
low T as expected for Fermi liquids, a nearly linear in- 
crease with increasing T for T < Tq, a maximum of R{T) 
exactly at To, and an R{T) decreasing with increasing T 
(and thus behaving similarly as in the case of incoherent 
scattering from magnetic impurities) for T > Tq. 



v„=o.3 T„ = 0.02 

— V„=0.4 T„ = 0.05 

— V„=0.5 T„ = 0.09 
V„=0.6 T„ = 0.16 

=0.7 T„ = 0.24 

=0.8 T„ = 0.30 

=0.9 T„ = 0.37 




FIG. 5. The resistivity of the periodic Anderson model is 
shown for various values of the hybridization parameter as a 
function of temperature measured in units of the band-width 
W. The total number of electrons per site is ntot = 2.2 and the 
number of / holes is n^"'" = 0.9. The value of the Coulomb 
correlation is U/W — 1 and the corresponding low-energy 
scale To is indicated in the figure. 
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As an example, we consider the model with more than 
one / electron per site. This would correspond to Yb 
systems in which the number of / holes is restricted by 
large Coulomb correlation to nj"'*^ < 1. The numeri- 
cal results are shown in Fig. where R{T) is plotted 
versus temperature for the model with the total number 
of electrons per site ritot = 2.2, the occupancy of the / 
holes n^"'-'^ — 0.9, and for different values of hybridiza- 
tion V. For fixed other parameters the hybridization V 
determines the low-temperature scale Tq, which is indi- 
cated in the figure. The calculations are performed for 
U/W = 1. Remarkably, as long as there is a maximum in 
i?(r), it is very close to Tq as determined by the / DOS 
criterion described above. For too large V and Tq (and 
the corresponding less strong mass enhancement) there 
is no longer a true maximum in R(T) but only a plateau 
behavior. 



-110 
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FIG. 6. The thermopower of the periodic Anderson model 
is shown for various values of the hybridization parameter as a 
function of temperature, measured in units of the band-width 
W. The total number of electrons per site is ritot = 2.2 and the 
number of / holes is rij"'" = 0.9. The value of the Coulomb 
correlation is U/W — 1 and the corresponding low-energy 
scale To is indicated in the figure. 



Corresponding results for the thermopower are shown 
in Fig. ^. Note that one obtains S{T) in its natural units 
(using ks/e « SGnV/K) from (16), as the "arbitrary" 
units (prefactor of (13)) cancels due to the quotient. The 
initial sign of S{T) depends on the slope of L{E) at the 
Fermi level and the temperature variation of S{T) re- 
flects the structure of the DOS within the Fermi win- 
dow. Fig. 1^ shows typical results for S{T) of Yb systems, 
with a large negative low-temperature peak. Obviously 
S{T) is strongly temperature-dependent and its absolute 
value is very large, namely of the magnitude 50 - 100 
/iV/K. It has an extremum (negative minimum in the 
plot) at a temperature Ti, which scales linearly with Tq, 
i.e., Ti = aTo. 

For the parameters used for the figure we have a w 0.5, 



but the exact value of a depends on other parameters (U, 
rif etc.) as well. This extremum in S(T) also exists in the 
situation when the resistivity R{T) exhibits no maximum 
but only a plateau behavior. The absolute value |S'(ri)| 
at the extremum can even increase with increasing Tq, 
i.e., according to this result it is not necessarily the most 
"heavy" fermion systems that exhibit the largest values 
of the TEP. In any case, in the low-temperature regime, 
T < Tq, this approach yields a TEP of the correct ab- 
solute magnitude and qualitative behavior, in particular 
an extremum, which is characteristic for the TEP exper- 
imentally observed in many heavy fermion systems. Of 
course for intermediate and high T the features due to 
the CF splitting cannot be reproduced within this SOPT 
treatment of the PAM, as only the twofold degenerate 
PAM was used, which has no higher CF-split / levels in- 
cluded. 



Single-impurity case 

The transport properties of dilute Ce and Yb alloys 
with the doubly degenerate CF ground state are obtained 
at temperatures below Tk from the single-impurity spin- 
1/2 Anderson model. We assume that the number of / 
electrons (holes) is slightly above (below) one for each Ce 
(Yb) impurity, and we consider an asymmetric model. 
The effective parameters of the SIAM are adjusted in 
such a way that the width of its Kondo resonance coin- 
cides with the Kondo temperature Tk of the CSM with 
CF splitting. 

The TEP of a spin-degenerate SIAM, in the absence 
of the non-resonant scattering channels, has been calcu- 
lated by various methods and is well understood. We 
are interested in the low-temperature behavior obtained 
for U/ttT ^ 1, where U denotes the f-f correlation and 
r = TTpF is the half-width of the /level due to the hy- 
bridization with conduction states. (Here, pp is the den- 
sity of conduction states at the Fermi level.) The numer- 
ical renormalizatiou-Ecpup calculations(NRG) of Costi, 
Hewson and ZlaticCJO (in what follows referred to as 
CHZ) and the perturbation theorycJ show that S{T) is 
closely related to the Kondo resonance. In the absence of 
any potential scattering the behavior of the TEP follows 
simply the temperature dependence of the it f-electron 
spectral function within the Fermi window. The results 
for the 4f^ electron (and analogousjeailts for the 4f^^ 
hole) can be summarized as foUowsJlS^LJ Close to T = 
the system is in the Fermi liquid (FL) regime, which is 
characterized by an asymmetric Kondo resonance of the 
width Tpc, centered above the Fermi energy, Ep. This 
leads to a typical Fermi liquid power law, 

SiT)^^^cot MEp), (17) 
3|e| Tk 

where •qo{Ep) is the resonant phase shift due to the scat- 
tering of conduction electrons on the / state. It is related 
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to Uf by the Friedel sum rule, rjo^Ep) = 7rn//2, so that 
for n/<l(?T,/>l) the initial slope of the TEP is pos- 
itive (negative). In heavy fermions, where temperatures 
of the order of Tk are easily accessible, the TEP grows 
rapidly with temperature (Tr^fcs/Slel — 284:fiV/K) and 
can assume giant values. At elevated temperatures the 
Fermi liquid behavior breaks down, because the spectral 
weight is transferred out of the low-energy region and the 
Kondo resonance disappears. For T > Tk, the system 
enters the loca- moment (LM) regime. Here, for n/ > 1, 
the maximum of the spectral function shifts from above 
to below Ep and the TEP is negative. For Uf < 1, we 
find that in the LM regime the TEP is positive. 

These simple features are modified in the presence of 
non-resonant scattering at the impurity site, because the 
interference between the resonant and the non-resonant 
channels leads to vertex corrections, which have drastic 
effects on the thermopower. These effects have been an- 
alyzed in detail in CHZ, and here we use that theory to 
study the changes in S{T) induced by a small variation 
in n^. Such changes in could be induced by pressure 
or doping (chemical pressure). The TEP is still given 
by the expression ( p^ ) , but the |eyaluation of the Kubo 
formula gives for L(ti^) the resulted 



Liu;) LoiiJ,T) 



ReG(w) . „ 
cos 2r]i - - — — — sm 2ryi 
Im G(uj) 



(18) 



Here, r]i is the non-resonant phase shift, pn — sm^(?7i) 
is the residual resistance due to the non-resonant scat- 
tering, and Lq{uj) is the transport relaxation time due to 
the resonant scattering. 



1 F 

—— - = Im G(tj), 

Lo(uj) IT 



(19) 



(for details of derivation see CHZ). The resonant Green's 
functionj-Ln Eq. ( pj| ) has been calculated by the NRG 
methodsLJ and the results were used to study the ther- 
mopower at the crossover from above to below Tk- How- 
ever, in the presence of CF splitting we do not need 
the resonant Green's function of the spin- 1/2 problem 
for temperatures much above Tk- In the Fermi liquid 
regime it-is sufficient to evaluate G{uj) by an approximate 
methodEa that interpolates between the 2nd-order weak- 
coupling expression (in what follows referred to as MR) 
and the exact atomic limit. This simple approximation 
is quite accurate for T < Tk', it gives the spectral func- 
tion and the thermoelectric power similar to the NRG 
results of CHZ even for U/ttT ^ 1. The approximation 
becomes unphysical for T 3> Tk but here we are not in- 
terested in such high temperatures where the spin-1/2 
model becomes inadequate anyway. The approximation 
is based on the Dyson equation, 



Gito) 



ei + iT — Uuf — T,{lo) ' 



(20) 



with ei determined iteratively in such a way that the 
number of electrons per spin. 



dw/(cj) ( -ilm G 



(21) 



equals Uf. The self-energy S(w) is defined by the Martin- 
Rodero interpolation, 



S(^) 



S(2)(t^) 



1- (i^Ii^)Ci±£i_£o^(2)(^))' 



(22) 



where, eg is the energy of the virtual bound state of the 
U ^ Anderson model with nf electrons, and I](2)(z) 
is the 2nd-order self-energy calculated with the unper- 
turbed propagators of this auxilliary U = model (for 
details see Ref. MR). 

The thermopower obtained by this approximation is 
plotted in Fig. ^ as a function of temperature and for 
several values of n/. We consider the case of strong corre- 
lation, U /ttF = 8, and assume a small non-resonant scat- 
tering, 771 = —0.1 (the effect of p„ on S{T) is neglected). 
The value of Tk (in units of ttF) is estimated from the / 
spectral function that we obtain from the Dyson equation 
(P0|). We use two different estimates: (i) we associate Tk 
with the full width at half the maximum (FWHM) of the 
Kondo resonance at T = 0, and (ii) we identify Tk as 
the half-width of the temperature-dependent part of the 
it f-electron spectral function at the chemical potential. 
The value of Tk obtained by either method is almost the 
same, and we find that the reduction of Uf enhances Tk- 
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FIG. 7. The thermopower of the single-impurity Anderson 
model is shown for various / occupation as a function of tem- 
perature, measured in units of Tk- The value of the Coulomb 
correlation is U/nT = 8 and the corresponding Kondo scale 
Tk is indicated in the figure. 



The Uf = \ curve corresponds to an extreme Kondo 
limit, where S{T) goes rapidly to a large negative peak 
at about T ~ Tk/^, and then decreases gradually but 
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remains negative within the physically relevant tempera- 
ture range T < Tk- The effect of the non- resonant phase 
shifts is very important, because at half- filling (n/ = 1) 
the model is electron-hole symmetric and in the absence 
of the non-resonant phase shifts the thermopower van- 
ishes. The n/ — 0.98 and rt/ — 0.99 curves are still 
Kondo-like, but the reduction of 770 from 7r/2 gives rise 
to a positive initial slope of S{T), such that S{T) de- 
velops a well defined positive maximum before decreas- 
ing sharply to negative values for temperatures between 
Tk/2 and Tk- Finally, the n/ = 0.96 and nj = 0.92 
curves remain positive for T < Tk- Further reduction of 
Uf does not bring any new qualitative features, except 
that the maximum of S{T) becomes higher and shifts to 
higher temperatures. By changing the other parameters, 
like the correlation or the non-resonant phase shifts, we 
find a similar set of curves, that could be classified in the 
same way as the curves shown in Fig. (0). 

The low-temperature thermopower results shown in 
Fig. (|^) exhibit all the features found in dilute Ce al- 
loys for T < Tx- Reflecting the curves on the 5* = axis, 
gives S{T) that correspond to the Yb alloys, with one 
magnetic / hole. 

Discussion of different cases (a) to (e) and conclusion 

The overall temperature dependence of the TEP of 
Ce and Yb intermetallics is found by interpolating be- 
tween the results for the Anderson model, valid below 
Tk, and the CSM results, valid above Tk- The essential 
features of S{T) are due to the Kondo effect modified 
by CF splitting, and various shapes are the consequence 
of different energy scales that characterize the low- and 
high-temperature behavior. The details, however, de- 
pend not just on the single parameter Tr-, but on other 
parameters of the model as well. Here we discuss, for 
simplicity, the thermopower due to the cerium ions with 
a single electron in the 4f^ configuration. The results per- 
taining to ytterbium ions with a single hole in the 4f^^ 
configuration are obtained by reflecting the S{T) on the 
temperature axis. 

In (a)- and (b)-type systems a large CF splitting of 
Ce / states generates two distinct energy scales, which 
are seen as two Kondo temperatures of the CS model, 
Tk and T^ ^ Tk- The high-temperature Kondo scale 
describes a sixfold degenerate CF multiplet with a large 
thermopower. Typically, T^ is between 30 and 100 K 
and the maximum value of the thermopower can exceed 
150 /iV/K. The low-temperature scale, deflned by Tr-, 
is of the order of 1 to 10 K. A twofold or fourfold de- 
generate CF ground state gives rise, below Tr-, to an 
additional structure in S{T)- In periodic systems, the 
onset of coherence reverses the sign of the thermopower 
below Tk and gives rise to a large negative peak around 
To/2 ~ Tk/I- The behavior of S{T) in the T = limit 
depends on the band flUing and the value of the f-d hy- 



bridization: for systems close to the electron-hole symme- 
try S{T) is positive for Uf <1 and negative for n/ > 1. 
Combining these results with the T < Tk results one ob- 
tains the thermopower of the (a) or (b) type. In dilute 
Ce alloys the sign of the low-temperature thermopower 
is determined by the non-resonant phase shifts. Type (b) 
behavior is obtained here for a ground CF-level that gives 
rise to a positive thermopower around Tx/2 and a nega- 
tive one above Tk- For sufficiently different values of Tk 
and T^ the S{T) has two well resolved positive peaks 
separated by a negative minimum. Changing the sign 
of the non-resonant phase shifts one reverses the sign of 
the low-temperature thermopower and finds the (a)-type 
thermopower: there is a deep negative minimum around 
Tk/2, followed by a positive maximum around T^ - The 
type- (a) and type-(b) systems can be classified as Kondo 
systems with large CF splitting. The crossover from the 
high- to the low-temperature fixed point is accompanied 
by the sign change of the TEP. 

The type-(c) behavior can arise in two ways. In the 
first case, the low- and high-temperature maxima over- 
lap, giving a two-hump structure. That is, the CF split- 
ting is not large enough to generate sufficiently differ- 
ent Kondo scales for the negative minimum to occur. In 
the second case, the CF splitting might be large but the 
potential scattering is also large enough to prevent the 
thermopower due to the CF ground state from changing 
sign for T < Tk- Thus, the crossover from the low- to 
high-temperature fixed point occurs at Tk without the 
change of sign of S{T). The two LM regimes are still well 
resolved but S{T) is always positive. 

The type-(d) thermopower is obtained for Kondo sys- 
tems with small CF splitting, such that the low- and 
high-temperature thermopower peaks overlap, and we 
see but a single peak with, perhaps, a shoulder on the 
low-temperature side. In the absence of the CF splitting 
there is just a single thermopower peak at about Tk/^,- 
The sixfold degenerate 4f^ electron of Ce ions has a large 
Kondo scale and huge positive thermopower. (The eight- 
fold degenerate 4f^'^ multiplet of Yb ions could lead to 
even larger negative thermopower.) 

In valence fluctuators, the / level is close to the chem- 
ical potential and Ce ions fluctuate between two con- 
figurations. The thermopower of such a system increases 
monotonically to a large positive (or negative) value. The 
characteristic energy scale is defined by the width of the 
/level, and S{T) might reach very large values. 

Finally, we remark that in Kondo systems with Ce ions 
the application of pressure (or positive chemical pressure) 
enhances Tk and transforms an (a)- or (b)-type ther- 
mopower to a (c)- or (d)-type one. The negative (chem- 
ical) pressure reduces Tk and gives rise to a reversed 
behavior. In Kondo systems with Yb ions the pressure 
and chemical pressure have the opposite effect than in Ce 
systems. We hope this simple explanation of the overall 
behavior of the thermopower of Ce and Yb systems to 
facilitate the search for optimal thermoelectrics. 
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